Model for pairing phase transition in atomic nuclei 
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I. INTRODUCTION 



Phase transitions in small systems are an important 
research topic. One of the most interesting problems 
in this context is probably the question of the existence 
and classification of a possible phase transition from a 
hadronic phase to a quark-gluon plasma in high energy 
physics. The answer to this question has far-reaching 
consequences into many other fields of research like, e.g., 
cosmology, since it has been argued that hadronization 
of the quark-gluon plasma should be a first-order phase 
transition in order to allow for possible supercooling 
and consequently the emergence of large scale inhomo- 
geneities in the cosmos within the inflationary big-bang 
model |]. 

In nuclear physics, two different phase transitions have 
been discussed in the literature. A first-order phase tran- 
sition has been reported in the multifragmentation of 
nuclei ||^, thought to be the analogous phenomenon in 
a finite system to a liquid-gas phase transition in the 
thermodynamical limit. A pivotal role in these studies 
is played by the presence of a convex intruder in the 
microcanonical entropy curve [^J^ . This leads to a neg- 
ative branch of the microcanonical heat capacity which 
is used as an indicator of a first-order phase transitions 
in small systems. Negative heat capacities have indeed 
been observed in the multifragmentation of atomic nu- 
clei, though the heat capacity curve has not been derived 
directly from the caloric curve, but by means of energy 
fluctuations . Another finding of a negative branch 
of the heat capacity curve has been in sodium clusters 
of 147 atoms [Q, indicating a possible first-order phase 
transition. On the other hand, it is not clear whether 
the observed negative heat capacities are simply due to 
the changing volume of the system under study that is 
progressively evaporating particles |Q. It has also been 



suggested that a negative heat capacity can be observed 
when studying a system in a metastable state |^ . In gen- 
eral, great care should be taken in the proper extraction 
of temperatures and other thermodynamical quantities 
of a multifragmenting system P,p^. 

A further difficulty in characterizing the order of a 
phase transition in a finite system using negative heat 
capacities arises from the fact that there are no coun- 
terparts in the thermodynamical limit. For this reason, 
no simple connection with the Ehrenfest classification of 
phase transitions can be made. Therefore, it is not clear 
that a phenomenon involving negative heat capacities in 
small systems is analogous to a first-order phase transi- 
tion in the thermodynamical limit. Another way to clas- 
sify phase transitions for finite systems has been proposed 
within the canonical ensemble. This classification scheme 
is based on the distribution of zeros (DOZ) of the canon- 
ical partition function in the complex temperature plane 
|]lT| and has successfully been applied to a model of mul- 
tifragmentation ||T^. This classification scheme reduces 
to the Ehrenfest classification in the thermodynamical 
limit and is applicable for first, second, and higher order 
phase transitions. 

The second phase transition discussed for atomic nu- 
clei has been anticipated for the transition from a phase 
with strong pairing correlations to a phase with weak 
pairing correlations Early schematic calculations 

have shown that pairing correlations can be quenched by 
temperature as well as by the Coriolis force in rapidly 
rotating nuclei |p^-|l7||. This makes the quenching of 
pairing correlations in atomic nuclei very similar to the 
breakdown of superfluidity in '^He (due to rapid rota- 
tion and/or temperature) or of superconductivity (due 
to external magnetic fields and/or temperature). Re- 
cently, structures in the heat capacity curve related to 
the quenching of pairing correlations have been obtained 
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within the relativistic mean field theory ||l^,|i9[ , the finite- 
temperature random phase approximation (RPA) 
the finite-temperature Hartree-Fock-Bogoliubov theory 
|2l|, and the shell-model Monte Carlo (SMMC) approach 
1 22-p6t. An S-shaped structure in the heat capacity curve 
could also be observed experimentally and has been 
interpreted as a fingerprint of a second-order phase tran- 
sition from a phase with strong pairing correlations to a 
phase with weak pairing correlations. Indeed, the anal- 
ogy of the quenching of pairing correlations in atomic 
nuclei with the breakdown of superfluidity in '^He and 
the breakdown of superconductivity suggests a second 
order phase transition and a schematic calculation might 
support this assumption . Interestingly, similar struc- 
tures of the heat capacity curve as observed for atomic 
nuclei in pTf have been seen in small metallic grains un- 
dergoing a second-order phase transition from a super- 
conductive to a normal conductive phase |^^, thereby 
supporting the analogous findings for atomic nuclei. On 
the other hand, breaking of nucleon pairs has been ex- 
perimentally shown to cause a series of convex intrud- 
ers in the microcanonical entropy curve of deformed rare 
earth nuclei leading to several negative branches 

of the microcanonical heat capacity. This finding might, 
in analogy to the discussion of nuclear multifragmenta- 
tion, be taken as an indicator of several first-order phase 
transitions. Interpreting the quenching of pairing corre- 
lations as a series of first-order phase transitions seems, 
however, physically unattractive and might rather sug- 
gest that a negative branch of the microcanonical heat 
capacity curve does not necessarily indicate the analo- 
gous phenomenon in a small system to a first-order phase 
transition in the thermodynamical limit. A similar con- 
clusion has been drawn by the authors of Ref. |^ . 

We will in this work develop a model for the atomic nu- 
cleus which allows us to investigate the occurrence and 
classification of the pairing phase transition in atomic 
nuclei as a function of particle number and other param- 
eters within the model. In order to do so, we will apply 
the classification scheme of Borrmann et al. 12 1 to the 
pairing phase transition. 



II. STATISTICAL ENSEMBLES 



The proposed nuclear model is formulated within the 
canonical ensemble theory although the microcanonical 
ensemble theory might, at a first glance, be the more ap- 
propriate ensemble for describing the nucleus below the 
particle threshold. At low excitation energies, the nucleus 
can be regarded as a closed system with respect to energy 
and particle exchange. However, the microcanonical en- 
semble theory, just because it requires a closed system, 
has an inherent conceptual problem when it comes to 
small systems of few particles. Within the microcanoni- 
cal ensemble theory one has 



(1) 



For macroscopic systems, it can be proven experimentally 
that T from this equation equals the thermodynamical 
temperature defined by primary and secondary temper- 
ature standards, like the triple point of water. Unfor- 
tunately, for small systems, this experimental proof can 
in principle not be performed due to the following dif- 
ficulty. Temperature measurements involve necessarily 
energy exchange between the system and the thermome- 
ter in order to bring the thermometer to the same tem- 
perature as the system under study, thus undermining 
the requirement of a closed system for a microcanoni- 
cal description. For macroscopic systems, this problem 
is solved by introducing a small thermometer compared 
to the system. For small systems of a few particles, in 
general, every thermometer will, due to its size, act as 
a heat bath when brought into contact with the system 
under study. Thus, temperature measurements on small, 
truly microcanonical systems seem virtually impossible 
and, therefore, it cannot be proven that the quantity T 
from Eq. corresponds to the thermodynamical tem- 
perature in these cases, although S and E are still well 
defined and can be measured. It should also be noted 
that all means of temperature measurements which are 
based on the application of the statistical theory, like mo- 
mentum distribution measurements, are not valid either, 
since the connection between the microcanonical statis- 
tical theory and thermodynamics is always based on Eq. 

The treatment of the nucleus in the canonical en- 
semble theory has also the advantage that a classifica- 
tion scheme for phase transitions exists which reduces 
to the Ehrenfest definition in the thermodynamical limit 
Hence, the difficulties mentioned in the Intro- 
duction can be overcome. Although it was suggested 
in [^Ql that the canonical ensemble theory should not 
be applied to first-order phase transitions, it has been 
shown that it gives convincing results in the case of nu- 
clear multifragmentation which is thought to be a 
first-order phase transition. It has also been objected 
that the Laplace transformation involved in the canoni- 
cal ensemble theory smears out the thermodynamical in- 
formation of the system under study and that the canon- 
ical partition function therefore contains less information 
than the microcanonical partition function We want 
to emphasize, however, that the partition function of the 
canonical ensemble, when known on all of the right-hand 
half of the complex temperature plane, certainly does 
carry as much thermodynamical information of the sys- 
tem as the microcanonical partition function. The reason 
for this is that the inverse Laplace transformation (from 
the canonical to the microcanonical partition function) is 
mathematically defined in a unique fashion and hence the 
microcanonical partition function can in principle be re- 
covered unanimously. However, since the inverse Laplace 
transformation requires integration on the complex tem- 
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perature plane along a parallel to the imaginary axis, one 
has to expect that knowledge of the canonical partition 
function over all of the right hand half of the complex 
temperature plane is pivotal for revealing the complete 
thermodynamical information of the system under study. 
The classification scheme which we intend to ap- 

ply in this work, is based on the DOZ of the canonical 
partition function and therefore complies with this re- 
quirement. 



III. MODEL 

The model in this work is a further development of 
the previously proposed model in ||3^,Q. We will first 
recapitulate those features of the old model which we 
have adopted in the present work and then describe the 
changes and additions. The basic idea behind our model 
is the assumption of a reservoir of nucleon pairs. These 
nucleon pairs can be broken and the unpaired nucleons 
are then promoted into an infinite, equidistant, doubly 
degenerated single-particle level scheme. The nucleon 
pairs in the reservoir do not interact with each other and 
are thought to occupy an infinitely degenerated ground 
state. The nucleons in the single particle level scheme 
do not interact with each other either, but they have to 
obey the Pauli principle. The model is depicted in Fig. 
|l|. The essential parameters of the model are: the num- 
ber of pairs in the reservoir at zero temperature N, the 
spacing of the single-particle level scheme e, and the en- 
ergy necessary to break a nucleon pair 2 A. Quenching 
of pairing correlations is introduced into this model by 
reducing the required energy to break a nucleon pair in 
the presence of unpaired nucleons. We assume that for 
every already broken nucleon pair, the energy to break a 
further nucleon pair is reduced by a factor r < 1. Thus, 
in the presence of p broken nucleon pairs, the energy to 
break the (p -I- l)th nucleon pair is reduced to 2 r^* A. We 
assume further that the single unpaired nucleon in an 
odd nucleus reduces the necessary energy for breaking 
the first nucleon pair of the same species to 2 v^A. 

The proposed model is supposed to display three dis- 
tinct phases: a phase where pairing correlations domi- 
nate, a phase where the pairs have essentially been bro- 
ken up and the system behaves like a non-interacting 
Fermi-gas dominated by Pauli-blocking, and a phase 
where the nucleons become diluted in a phase space 
much larger in size (in terms of the number of realiza- 
tions of thermal excitations) than the number of nucleons 
present. In this third phase, the Pauli blocking can be ne- 
glected and the system resembles mostly a classical ideal 
gas. We will, in the following, denote the three phases 



as the paired phase, the unpaired phase, and the quasi- 
classical phase. The three distinct phases of the system 
makes it possible to investigate two different phase tran- 
sitions with our model. 

It has been shown previously that the partition 
function for n unpaired nucleons in an infinite, equidis- 
tant, doubly degenerated single-particle level scheme can 
be written as 



(n-l)/2 

E 

(n-2)/2 

E 

i=0 



2 ^n-i 



= odd 



2 4-. ^1 + (4/2)' 



(2) 
(3) 



where the zj are the partition functions for i spin-up nu- 
cleons in an infinite, equidistant, non-degenerated single- 
particle level scheme. The zj can be expressed by means 
of the recursive formula |32 



exp(— ze/3) 



l-exp(-(i-hl)e/3)' 



(4) 
(5) 



where /3 = According to our model, in order to 

promote n nucleons from the reservoir into the single- 
particle level scheme, one has to break p = n/2 nucleon 
pairs which requires the energy 



2A + 2rA 



2r^A 



2rP-^A = 2A 



1 



1 - r 



(6) 



in the even case. In the odd case, one has to multiply this 
energy by ^/r in order to take into account the quenching 
effect of the odd nucleon. Hence, the partition function 
for n unpaired nucleons in our model becomes 



Zn = z„ exp(-2 ^/rA 
Zn = Zn exp(-2A-"'^ 



1 



,(n-l)/2 



1 - 

j,n/2 



1 



r 



even 



n ~ odd (7) 
(8) 



In order to obtain the partition function for the complete 
system, one has to sum over all possible numbers of un- 
paired nucleons in the single-particle level scheme. Thus, 
the partition function of our model is finally 



Z„ 



N+l 

Z2n-1 



71=1 

N 



z 



pair 



- Z2„ 



odd 



even. 



(9) 
(10) 



ri=0 



^Here, and in the following, we set Boltzmann's constant 
cb = 1. 
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At this point, the average number of unpaired particles 
can be calculated in a straightforward manner 



N+l 



N+1 



n=l 

N 

''Z2n/'^Z2n evcn. 



odd 



n=l 
N 

n=0 



(11) 

(12) 



n=0 



The number of unpaired particles will fluctuate, and the 
fluctuations can be calculated by 



5n = \/ {n?) — (n)2, 
where {it?) is given by 

N+l N+l 



(13) 



n=l 



= Z2n-l/ ^2n-l odd (14) 

(15) 



ri=l 

N N 
n=0 ri=0 



) = Z^n Z2nl y . ^2n even 



Since there are protons and neutrons in the nucleus, 
one has to create two partition functions Zpair, one for 
every nucleon species, and multiply them with each other. 
Here, one can introduce different parameters e. A, r, and 
N for protons and neutrons. For the sake of simplicity, 
we will, however, in this work not make use of this pos- 
sibility^. In addition to the model partition function, we 
have in the previous work | |3^ ] also assumed a rotational 
and vibrational partition function 



12 

E 

7=0,2,4.. 
1 



exp(-Arig/(/ + 1)13) 



i/=0 



(16) 



(17) 



with Ai-ig being the rigid-body rotational parameter and 
Wvib the vibrational frequency of the nucleus under study. 
The nuclear partition function becomes therefore in the 
most general case 



^ ^pair ^pair ^lot ^vib 



(18) 



where we will use in this work equal proton and neutron 
partition functions Z^^^^ and .^pair- We would like to 
stress that for r = 1, the present model is approximately 
equal to the previous model js^js^, whereas the possi- 
bility of r < 1, i.e., the quenching of pairing correlations, 
is a new feature of the present model. 



Although the model is purely phenomenological, we 
would like to emphasize that it can be closely related to 
a microscopical model. If one assumes a Hamiltonian of 
the form 



H — Hjr + Hi, + Hrot + -ffvib 



(19) 



with 



Hn,u = eY<^ici. - \G\ E 4aLa-AaA, (20) 

K K,A>0 

i.e., a simple single-particle plus pairing Hamiltonian for 
protons and neutrons, and conventional rotational and 
vibrational Hamiltonians i?rot and H^ih for the collective 
modes, similar results as in this work can be achieved 
1^,0. The advantage of the present, phenomenological 
model over the microscopic model is that, in general, a 
much larger number of particles and single-particle levels 
can be taken into account. 



IV. MODEL PROPERTIES 

It has already been shown that the previous model 
p2| , ^ describes well the level density and heat capacity 
in the ^^^Dy nucleus. Further, it could reproduce 3^ ] 
the anchor points of the level density curve proposed in 
p2[ for several other mid-shell nuclei. However, in order 
to do so, one had to assume values for e which were off 
general estimates of the shell model. Also, it was not 
possible to reproduce the theoretical heat capacity curve 
of iron nuclei as calculated in ||2^ within the SMMC ap- 
proach. Since in the new model an additional parameter 
is introduced, we have the hope to use more reasonable 
values of e and still be able to describe the experimental 
data with good accuracy. 

The parameters for the present nuclear model are taken 
from the following systematics [B5| 



37r2 
A 



MeV, 



(21) 



where the level density parameter a = Tr^/Se is related to 
the mass number Ahy a — A/9MeV~^ . The pairing-gap 
parameter A can also be related to the mass number A 
by 



12 

A = ^ MeV. 

Va 



(22) 



The rigid-body rotational parameter Aiig can be ex- 
pressed in terms of the rigid moment of inertia Orig by 



^This means that the even-odd and the odd-even systems 
will be equal. We will, therefore, in this work use the term 
'odd system' to denote either of the two. 
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(23) 



piiE)) 



Assuming the nucleus to be a rotational symmetric ellip- 
soid of constant density, the rigid moment of inertia can 
be expressed by the nuclear mass M and mean radius R 
as Grig = 2/bME?. The nuclear mean radius can then 
again be related to the mass number by i? = 1.24 fm A^/'^. 
The vibrational frequency of the nucleus can be taken di- 
rectly from spectroscopic data [Q. For the parameter r 
of our model, which governs the evolution of the pairing 
gap A with temperature, no systematics exist. Theoreti- 
cal studies suggest a quenching of the pairing energy per 
broken pair in the order of 0.56 in the mass 190 region | 3^ ] 
yielding a A(r) curve comparable to results from mod- 
ern RPA calculations j38| (see Fig. Probably r will 
increase with mass number. This effect is, however, not 
taken into account in this work and we adopt r — 0.56 for 
all nuclei studied. Finally, for the last parameter N of the 
model, we assume ten neutron and ten proton pairs in the 
reservoir at zero temperature in most of this work. This 
is an increase over the previous calculations [ ^2|j3^ where 
a total number of ten nucleon pairs was assumed. The 
increase of the number of nucleon pairs in the reservoir 
allows us to map thermodynamical quantities of nuclei 
up to higher temperatures of, e.g., 1.5 MeV in dyspro- 
sium. A summary of the applied parameters is listed in 
Table |. 

In order to compare the model to experimental data, 
one has to calculate experimentally observed quantities 
like the free energy F , the entropy S, the caloric curve 
{E{T)), the heat capacity Cy, and the nuclear level den- 
sity p{{E)) from the canonical partition function. While 
the first four quantities can be calculated in a straight- 
forward fashion by the following expressions 



exp(S') 



(28) 



F 
S 



-TlnZ 
'dF 
df 



(E) ^F + ST 



(24) 
(25) 
(26) 
(27) 



the nuclear level density should, in principle, be obtained 
by an inverse Laplace transformation of the canonical 
partition function. We will, however, in this work make 
use of the saddle-point approximation |23[ 



which gives satisfactory results for the nuclear level den- 
sity ||33| . Finally, we will calculate the average number of 
unpaired nucleons (n), the fluctuation of this number Sn, 
and the effective, temperature-dependent pairing energy 
A(T) =rW2)A(0). 

In Fig. ^, thermodynamical quantities derived from the 
model are shown using parameters for ^^^Dy. The quan- 
tities S, {n), Cv and A(T) are coinciding for the even- 
even, odd, and odd-odd systems above temperatures of 
0.8-1 MeV, indicating that the pairing correlations are 
already completely quenched at these temperatures. In- 
deed, the heat capacity curve displays a linear depen- 
dence on temperature between T = 1 and 1.5 MeV, which 
is the expected temperature dependence for a Fermi gas. 
The fluctuation of the number of unpaired particles Sn 
increases strongly around T ~ 0.5 MeV, which coincides 
with a strong decrease of the effective pairing energy 
A(T) and the S-shape of the heat capacity curve, first 
observed experimentally in The latter three obser- 
vations can be interpreted as signatures of a second-order 
phase-transition like phenomenon in nuclei. 

On the right-hand panels of Fig. |^, the quantities {n) 
and C'v can be seen to split up again for the even-even, 
odd, and odd-odd systems at temperatures above 1.5- 
2 MeV. At these temperatures, all nucleons start to be- 
come excited, the reservoir is rapidly being exhausted 
and the system can no longer absorb energy at the same 
rate without getting heated up rapidly. Therefore, the 
heat capacity decreases until it reaches the limit of 1 fee 
per particle'^ and the number of unpaired particles will 
increase at a progressively smaller pace to the limit of 40, 
41, and 42 for the even-even, odd, and odd-odd systems, 
respectively. This is the signal of a possible second phase 
transition inherent in the model from an unpaired phase 
to a quasi-classical phase. It is not expected that this 
second phase transition is present in the experimental 
data since the nuclei under study in this work have more 
particles than are included in the model space. There- 
fore, the second phase transition in the model merely 
reflects the fact that the model space is still too small in 
order to describe heavy nuclei adequately at the highest 
temperatures*. Some kind of signal of the exhaustion of 
the model space will, however, be present in all types of 



^This corresponds to the value for a two-dimensional ideal, 
classical gas. 

*In any case, atomic nuclei become unbound, undergo multi- 
fragmentation and therefore loose their identities at the criti- 
cal temperature for the liquid-gas phase transition. This will 
therefore be the temperature beyond which our model be- 
comes unrealistic, even when including all nucleons of the 
nucleus. 
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nuclear model calculations which do not include all nucle- 
ons or the complete single-pa rticle level scheme. We will 
therefore investigate in Sect. VB whether and how the 
truncation of the model space interferes with the signal 
of the pairing phase transition. 

Figure || shows how the model compares to the level- 
density anchor points of Rcf. and the measured level 
densities of two dysprosium isotopes |2^. The data 
points for the dysprosium isotopes have been modified 
slightly for the purpose of this work. While in the pre- 
vious works [ p7|pS| ] a back-shifted level-density formula 
with the von Egidy parameterization pO| ] was used for 
normalization and extrapolation to the neutron-binding 
energy and beyond, in this work, we will make use of a 
simple, shifted level-density formula 



P = f 



exp(2> 



12V0.1776ai/2[/3/2^i/3 



(29) 



with U ^E~A,a^A/9 MeV'^ and A = 12/ VA MeV 
for consistency. The parameterization for the spin cut-off 
parameter a is taken from Ref. [ p| and the normalization 
parameter / is chosen such that the level-density formula 
intersects the data point from neutron resonance spacing 
Q. It is satisfying that the model reproduces the data 
points well without any fine tuning of the parameters 
in the different mass regions off closed shells. Only in 
the case of the iron mass region, the single particle level 
spacing had to be increased from the systematic value of 
e = 517 keV to 800 keV in order to describe the data. The 
shortcomings of the model to reproduce the iron data us- 
ing parameters from the systematics is most likely due to 
the presence of the double shell closure in the neighbor- 
ing ^^Ni isotope. Irregularities of both, the single-particle 
level spacing (or more precisely the closely related level- 
density parameter a) and the pairing-gap parameter A in 
the vicinity of closed shells are a well-known fact . We 
would like to stress that a better fit to the experimental 
data might certainly be achieved by fine tuning of the pa- 
rameters for protons and neutrons for every mass region. 
However, the possibility of using systematic parameters 
for the model highlights its general applicability. 

In Fig. ^ the model predictions are compared to pub- 
lished canonical thermodynamical quantities. For the 
iron nuclei these data are recent SMMC calculations , 
while for the dysprosium nuclei, semi-experimental data 
exist p7| . The comparison is made for canonical heat- 
capacity curves, since these curves are effectively second 
derivatives of the canonical partition functions and there- 
fore more sensitive to structural changes in the nucleus 
than any other thermodynamical quantities. 

Although the steepness of the S-shapc differs somewhat 
and the absolute value at the plateau is off by ~ 30%, the 
overall agreement of the model calculation for the dyspro- 
sium nuclei with the semi-experimental data is satisfying. 
Especially the temperature where the steepness of the S- 
curve reaches the maximum, i.e., the critical temperature 
of the pairing phase transition, is well reproduced. For 



the iron nuclei, the agreement between our model and 
the SMMC calculations is not so good. Unfortunately, no 
experimental data exists to compare the two models to, 
therefore it is difficult to judge which is the more realistic 
one. While the present model is very schematic, it offers 
the possibility to include infinitely many single-particle 
levels and much more nucleons than the SMMC calcula- 
tions. Especially the necessity of a finite single-particle 
level scheme in SMMC calculations will cause the heat ca- 
pacity curve to approach zero for high temperatures and 
it is not clear at which temperature this effect takes over 
in the data of Ref. ||6|. On the other hand, the SMMC 
calculation can incorporate shell-effects in a more natu- 
ral way than by simply adjusting the two parameters e 
and A. Especially the use of a realistic nucleon-nucleon 
interaction and the correct single particle energies for the 
complete {pf + 099/2) shell might change the picture sig- 
nificantly at low temperatures compared to the present 
schematic model. 

In conclusion, it can be stated that the present model, 
using a systematic parameter set, reproduces fairly well 
the available experimental data, but a discrepancy with 
SMMC calculations in the iron-mass region exists. Since 
the agreement with experiment is good in the rare earth 
mass region, where also experimental data is abundant, 
we are confident that the model can be applied to classify 
the order of the pairing phase transition with good accu- 
racy in deformed rare earth nuclei if not in most mid-shell 
nuclei. 



V. CLASSIFICATION OF PHASE TRANSITIONS 
A. Classification scheme 

The classification scheme for phase transitions is based 
entirely on ||l^,|l^ and is only slightly modified for the 
purpose of this work. We will therefore only recapitulate 
briefly the main features of the classification scheme. The 
scheme relies on the DOZ and derives three quantities, 
where two of them classify the order of a possible phase 
transition and the third one reflects the discreteness of 
the system under study. 

First, we define the inverse complex temperature 



B = P + iT, 



(30) 



where f3 — 1/T as usual and r denotes the imaginary 
part of the inverse complex temperature which is mea- 
sured in MeV~^ in this work. The zeros of the canonical 
partition function typically line up on curves through the 
complex temperature plane. Whenever such a curve in- 
tersects the real axis, a phase transition in the system 
might occur. Figure ^ shows a typical example with four 
zeros on a curve approaching the real axis. The zeros are 
denoted by (/3j , Tj) with j = 1 . . . 4 and j increasing with 
increasing distance from the real axis. 
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The angle v, being the first quantity of interest, is cal- 
culated by 



V — arctan 



T2 - Tl 



(31) 



Further, the average inverse distance between zeros can 
be calculated as 



1 



1 



dj-i 



(32) 



where do 



'i+i 



r,)^. The function 



$ can then be approximated in the vicinity of the real 
axis by a power law of Tj only 



$(r,)(xr;. 



(33) 



such that the second quantity of interest, a, can be cal- 
culated by means of 



a = 



_ ln$(T3) - ln$(r2) 



In Ts — In T2 



(34) 



One can easily see that in order to calculate a, one has 
to determine at least the first four zeros closest to the 
real axis. With a slight redefinition of the function $ by 



(35) 



where Tj = (tj + rj_|_i)/2, only the first three zeros have 
to be known. The order of the phase transition is now 
completely determined by the two quantities v and a. 
When a<Oora = i/ = the system exhibits a first- 
order phase transition. For < a < 1 and z/ = or 

7^ 0, a second-order phase transition takes place. For 
a > 1 the phase transition is of higher order. 

The third parameter ti indicates the discreteness of 
the system. More precisely, the quantity Ti/fi equals the 
time after which an ensemble of equal systems loses its 
memory. Only if ti — > for increasing particle numbers, 
one has a corresponding phase transition in the Ehrenfest 
sense for the thermodynamical limit. For finite systems, 
one can then speak of a phase transition in the gener- 
alized Ehrenfest sense. Finally, (3i might be defined as 
the inverse of the critical temperature of the phase tran- 
sition. In the following, the classification scheme will be 
applied to the DOZ of mid-shell nuclei calculated within 
the proposed nuclear model. 



B. Results 

First, we investigate the effect of the numbers of pairs 
in the reservoir at zero temperature on the DOZ. On the 
upper left panel in Fig. ^ the calculation is done for one 
single pair in the reservoir. Certainly, no phase transition 
can be claimed, but the positions of the zeros already in- 
dicate at which temperatures the pairing phase transition 



will take place. Also, the effect of the Pauli blocking by 
the single unpaired nucleon is visible, shifting the zero of 
the odd system to higher temperatures. For two pairs in 
the reservoir, two zeros for the even and odd system can 
be seen. The zeros are positioned on slanted lines where 
the zeros of the odd system are again shifted to higher 
temperatures due to the Pauli blocking. The same is true 
for three pairs in the reservoir. According to the previ- 
ous subsection, three zeros in a line can already be used 
to classify a phase transition. In the present case, since 
a is negative, a first-order phase transition is indicated. 
Due to the small size of the system, however, one cannot 
attribute this phase transition simply to the quenching 
of pairing correlations. Rather, one has to acknowledge 
the fact that the pairing phase transition coincides with 
the phase transition from an unpaired to a quasi-classical 
system, giving rise to a disturbed signal for the order of 
the phase transition. This becomes clearer in the case 
of four pairs in the reservoir. Here, one might already 
discuss two separate phase transitions, since one can dis- 
tinguish two different curves of zeros approaching the real 
axis. While the curves to the right might be attributed 
to the pairing phase transition, the curves to the left are 
indicating the transition to a quasi-classical phase. Still, 
since the two curves are sharing zeros one cannot expect 
to get a clean signal for the classification of the two phase 
transitions and indeed for both curves one obtains neg- 
ative values of a, indicating two first order phase tran- 
sitions, a result which is at variance with expectations. 
First the calculation with seven pairs separates the two 
phase transitions enough in temperature to get clean sig- 
nals. Here, one finds negative values of a for the tran- 
sitions to a quasi-classical phase, indicating a first-order 
phase transition. Interestingly, this phase transition does 
show only a very small odd-even effect, which can in to- 
tal be attributed to the by one higher particle number 
in the odd case. The pairing phase transition, however, 
has a large odd-even dependence and yields a ~ 0.86 for 
the even case, indicating a second-order phase transition 
and a value greater than one for the odd case, indicat- 
ing a higher order phase transition. This result remains 
stable for larger particle numbers. The values of v are 
consistent with the conclusions drawn above. 

We will now investigate the evolution of the parameters 
/3i, i.e., the inverse critical temperatures of the respective 
phase transitions, as function of N , i.e., the size of the 
model space. In the upper left panel of Fig. 0, it is shown 
that Pi for the phase transition from the unpaired to the 
quasi-classical phase obeys a power law of the form aN^ 
with b ^ —1. This means that in the thermodynami- 
cal limit the critical temperature of this phase transition 
will approach infinity and thus, the phase transition can- 
not be observed in infinitely large systems. This is not 
surprising, since the very nature of this phase transition, 
which is a transition from a Fermi gas to a quasi-classical 
phase where the nucleons are diluted in phase space, obvi- 
ously requires increasingly higher temperatures the more 
nucleons are present. Here, we want to stress once more 
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that this phase transition is a true phase transition in- 
herent in the model and might even be reahzed in nature 
by some physical system. However, when applying the 
model to heavy nuclei which are much larger than the 
model space, Pi merely indicates at which (inverse) tem- 
perature the reservoir is being exhausted and that the 
model space is too small to describe heavy nuclei ade- 
quately below Pi. Eventually, the temperature beyond 
which our model becomes unrealistic, even when includ- 
ing all nucleons of the nucleus in the model space, is 
given by the critical temperature for the liquid-gas phase 
transition (see also 

On the lower left panel of Fig. the evolution of Pi for 
the pairing phase transition as function of the size of the 
model space is shown. The inverse critical temperature 
Pi of the pairing phase transition is independent of the 
size of the model space as it is expected, but it shows a 
pronounced odd-even effect due to the Pauli blocking of 
the single unpaired nucleon. 

Finally, we will investigate the parameters ri as func- 
tion of N in order to gain information on whether the 
phase transitions observed in the finite system under 
study correspond to true phase transitions in the ther- 
modynamical limit. Unfortunately, in the case of the 
first-order phase transition from the unpaired phase to 
the quasi-classical phase. Pi and therefore also ti show a 
trivial scaling with N which is clouding the relevant, ad- 
ditional N dependence of ri. Therefore, one has to inves- 
tigate the dimensionless parameter ti/Pi for a residual, 
non-trivial N dependence in order to be able to conclude 
on the reality of the phase transition in the generalized 
Ehrenfest sense for the finite system. On the upper right 
panel of Fig. |^ the dependence of Ti/Pi on A'' for the 
phase transition from the unpaired to the quasi-classical 
phase is shown to obey the power law aN'^ with b ^ —0.5. 
Thus, the dimensionless parameter ti/Pi is shown to ap- 
proach the real axis in the thermodynamical limit and the 
modeled phase transition satisfies the generalized Ehren- 
fest definition of a phase transition. 

For the pairing phase transition, one could in princi- 
ple investigate ti directly as function of N , since Pi is 
independent of N and thus, no trivial scaling of ti with 
N will alter the results. However, to treat everything on 
an equal footing, also in this case, we have investigated 
the dimensionless parameter ti/Pi. Since Ti/Pi shows 
no sign of approaching the real axis (see the lower right 
panel of Fig. |7|) , one has to ask the question whether the 
modeled pairing phase transition corresponds to any true 
phase transition in the thermodynamical limit. 

In order to answer this question, we introduce the no- 
tion that, in the case of the pairing phase transition, the 
size of the system is not determined by the size of the 
reservoir N, but rather by the mass number A and sub- 
sequently by the mass-number-dependent parameters e 
and A. Specifically, the parameter g = A/e scales as vC4 
and should thus be increased in order to model a larger 
system. This is in sharp contrast to the analysis in [ p8[ , 
where the authors assume that the simple increase of the 



number of particles N in the calculation should model 
heavier nuclei. 

Therefore, in order to investigate the evolution of ti 
for the pairing phase transition with respect to the rel- 
evant size of the system, i.e., the mass number A, we 
have calculated the DOZ for nuclei in four different mass 
regions (see Fig. Obviously, the pairing phase transi- 
tion can only be observed for nuclei with mass numbers 
A > 100 within the presented model, since for lighter 
nuclei, the curve of zeros relevant for the quenching of 
pairing correlations points away from the real axis, and 
therefore, one cannot claim, within the model, the exis- 
tence of a pairing phase transition for light nuclei. On the 
left-hand panel of Fig. ^, the inverse critical temperature 
of the pairing phase transition is shown to scale approx- 
imately with A"^^^ in the rare-earth mass region, which 
is expected when assuming Tc = A/ Si, where 5*1 is the 
single-particle entropy p2[ . The iron data do not quite 
follow this trend, since we have used a single-particle level 
spacing for the iron nuclei outside the global systematics 
of Eq. (H^), see also Table |. Again, by investigating the 
dimensionless parameter ti /Pi of the pairing phase tran- 
sition, the trivial mass dependence is taken out and the 
residual dependence of Ti/Pi on A will tell if the mod- 
eled pairing phase transition in the finite system under 
study is a phase transition in the generalized Ehrenfest 
sense. On the right-hand panel of Fig. ^ this residual 
A dependence of ri / Pi has been fitted by the power law 
aA'^ with b ^ —1/3, again omitting the iron data due to 
their unsystematic value of e. Thus, ti/Pi will approach 
the real axis for large mass numbers. Therefore, one can 
conclude that the modeled pairing phase transition is a 
true phase transition in the generalized Ehrenfest sense, 
provided that one investigates this phase transition with 
respect to the relevant size of the system, which is the 
mass number A rather than the size of the model space 
N. As a final remark, it is interesting to notice that the 
odd-even difference of Pi and ti /Pi becomes smaller with 
increasing mass number. Thus, in the thermodynamical 
limit, it will not matter if the system consists of an even 
or an odd number of particles, as it is expected. 



VI. CONCLUSIONS 

We have developed a model which allows the investiga- 
tion of phase transitions in a small system of two or more 
particles. The anticipated phase transition from a paired 
to an unpaired system in rare earth nuclei could be pro- 
duced. The order of the modeled pairing phase transition 
could be determined for even systems to be of second or- 
der and for odd systems to be of higher order. A second 
phase transition, which is of first order, is inherent in the 
model and connects the unpaired and a quasi-classical 
phase. The latter phase is characterized by the dilution 
of particles in phase space and might therefore be encoun- 
tered in systems where the phase space is much larger in 
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size (in terms of the number of realizations of thermal 
excitations) than the number of particles present. This 
second phase transition is a true phase transition and 
might be realized in nature by some physical system. In 
the case of modeling heavy nuclei with particle numbers 
much larger than the model space, however, it reflects 
merely the fact that the model space is still too small 
in order to describe those nuclei adequately at high tem- 
peratures (within the limits discussed in It has been 
shown that the pairing phase transition and the signal 
of the exhaustion of the finite model space have to be 
carefully separated in temperature in order to obtain a 
clear signature of the pairing phase transition in atomic 
nuclei. The presented model is able to take into account 
a sufficient number of single-particle levels and nucleons 
to allow a good separation and to display conclusive ev- 
idence of the pairing phase transition in heavy mid-shell 
nuclei. 
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TABLE I. The model parameters have been calculated by means of e = 30 MeY/A, A = 12 MeV/v^, and 
^rig = 34 MeV A~^/^ with the exception of e for ^*Fe (see text). All values are rounded to integer numbers. The values 



for fifjJvib have been 


taken from 


spectroscopic data of Ref. 


[36]. The value of r = 


0.56 is suggested for the mass 


190 region by 


the calculations of Ref. [37] and adopted for all four nuclei 


in this work. 






Nucleus 


e [keV] 


A [keV] 


Arig [keV] 


ftwvib [MeV] 


r 


^«Fe 


800 


1576 


39 


2.0 


0.56 


106p^ 


283 


1166 


14 


1.4 


0.56 




185 


943 


7 


0.9 


0.56 


2:U|. 


128 


781 


] 


0.8 


0.56 
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FIG. 1. Schematic representation of the nuclear model. For details, see text. 
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FIG. 2. Thormodynamical quantities calculated within the nuclear model. The solid, dashed, and dash-dotted lines are for 
the even-even, odd, and odd-odd system, respectively. For A(T), only the even and odd case for one nucleon species is shown. 
The long tail of the A(T) curves at high temperatures agrees well with modern RPA calculations [38]. The parameters are 
those for ^'^^Dy. 
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FIG. 3. Comparison of calculated and experimental level density. The open symbols are the anchor points of Ref. [32], the 
full symbols are the experimental level densities of Ref. [27] using a slightly different normalization (see text). 
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FIG. 4. Comparison of heat capacity curves of the present model and of other publications. The open symbols are scanned 
from Ref. [26], the data for the dysprosium nuclei are from Ref. [27] using a slightly different normalization (see text). 
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FIG. 5. Definition of important quantities for the classification of phase transitions in finite systems. 
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FIG. 6. DOZ for different numbers of pairs in the reservoir. The full symbols are for the even system, the open symbols for 
the odd system. The parameters of the model are those for '^'''^Dy. There exist zeros with T > 2 MeV"^ and/or /3 < MeV"\ 
but these are not of interest in this context. 



17 



> 

0) 




1.6 r- 



^ 1.55 

1 .5 
1 .45 

1 .4 
1 .35 

1 .3 
1 .25 

1 .2 
1.15 

1.1 



o 



10 



o 



20 



30 




0.74 
0.72 

0.7 
0.68 
0.66 
0.64 
0.62 

0.6 
0.58 
0.56 



40 



30 40 



o o 



• • • 



10 20 

number of particles 



30 40 



FIG. 7. Evolution of the parameters /3i and ri//9i as function of the size of the model space for the two phase transitions. 
Upper panels: the phase transition from the unpaired to the quasi-classical phase. Lower panels: the pairing phase transition. 
The full symbols in the lower panels are for the even system, the open symbols are for the odd system. The curves on the 
upper panels are fits to the data by a simple power law. The odd and even systems align on identical curves. 
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FIG. 8. DOZ for nuclei in four different mass regions. The full symbols arc for the oven systems, the open symbols stand 
for the odd systems. The line of zeros approaching the real axis perpendicularly corresponds to the exhaustion of the finite 
model space and is not relevant in the discussion. Only the more or less horizontal lines of zeros reflect the quenching of pairing 
correlations. 
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FIG. 9. Evolution of the parameters /3i and ti//3i of the pairing phase transition with respect to the mass number A, i.e., 
the relevant size of the system under study. The full symbols are for the even systems, the open symbols stand for the odd 
systems. The points for ''^'''''''Fe are not taken into account in the fit, since for the iron data, a single-particle level spacing e off 
the global systematics is used. 
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